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Abstract —In this paper we study the problem of identifying 
active suhcarriers in an OFDM signal from compressive measure¬ 
ments sampled at suh-Nyqulst rate. The problem is of importance 
in Cognitive Radio systems when secondary users (SUs) are 
looking for available spectrum opportunities to communicate over 
them while sensing at Nyquist rate sampling can be costly or 
even impractical in case of very wide bandwidth. We first study 
the effect of timing offset and derive the necessary and sufficient 
conditions for signal recovery in the oracle-assisted case when the 
true active sub-carriers are assumed known. Then we propose an 
Orthogonal Matching Pursuit (OMP)-based joint sparse recovery 
method for identifying active subcarriers when the timing offset 
is known. Finally we extend the problem to the case of unknown 
timing offset and develop a joint dictionary learning and sparse 
approximation algorithm, where in the dictionary learning phase 
the timing offset is estimated and in the sparse approximation 
phase active subcarriers are identified. The obtained results 
demonstrate that active subcarrier identification can be carried 
out reliably, by using the developed framework. 

1. Introduction 

Cognitive radio systems are emerging to respond to the 
ever-increasing demand for higher data rates and new wireless 
services by solving the spectrum underutilization problem [T], 
El. The most vital task in CR systems is Spectrum Sensing 
(SS) defined as identifying spectrum holes by sensing the radio 
spectrum and utilizing them without causing interference to 
primary users (PUs) 0. Of special interest in this regard 
is sensing of OFDM signals a, 0. OFDM is one of the 
most effective multicarrier techniques for broadband wireless 
communications which is employed by many of the current 
and emerging wireless technologies. 

On the other hand, due to the limitations of today’s 
analog-to-digital converter (ADC) circuits which cannot sup¬ 
port very high bandwidth and need excessive memory and 
prohibitive energy costs for implementing digital signal pro¬ 
cessing systems ||6l, it may be very costly and even impractical 
to sense the signal based on Nyquist-rate samples. This has 
motivated researchers to study sub-Nyquist methods for wide¬ 
band spectrum sensing in CR networks; see, e.g., m-M- 
However the number of works focusing on OFDM detection 
from sub-Nyquist samples is very limited. In lITSl . iflhl the 
authors proposed methods for detection of OFDM signals, 
but because both of these methods need a specific embedded 
signature inside the signal for the detection algorithm to work, 
they are suitable for detecting a cognitive network El rather 
than detecting PUs who are not willing to assist in signal 
detection process by embedding signatures in their signal. 
Besides, these methods as well as the one proposed in Cl 


do not address the problem of active subcarrier identihcation 
and rather focus on the detection of the OFDM signal as a 
whole. The works in Eol, El reconstruct the power spectral 
density (PSD) of a signal from its sub-Nyquist measurements 
which can be used for detecting used subcarriers of OFDM 
signal. The problem with these methods is that they are 
designed under the assumption of perfect timing information 
and therefore their performance degrades in the presence of 
unknown timing offset. Furthermore they need relatively many 
samples to compute the PSD of sub-Nyquist signal and then 
reconstruct the one of the original signal. 

Our contribution: In this paper, we propose a method for 
identihcation of active subcarriers of an OFDM signal with 
unknown timing offset. We assume that the OFDM signal is 
sampled by an Analog-to-Information converter (AIC) El, 
El at sub-Nyquist rates. We then model the system as a 
multiplication of a redundant dictionary which depends on the 
timing offset and a sparse vector whose non-empty entries 
correspond to the active subcarriers of OFDM signal. The 
effect of timing offset is studied and necessary and sufficient 
conditions for signal recovery in the oracle assisted case 
when the support of sparse vector is assumed known are 
provided. This is of great importance as each iteration of the 
proposed greedy algorithm involves a signal recovery step for 
the vector support estimated in that iteration. We then propose 
an algorithm that alternates between two steps: a dictionary 
learning step during which the timing offset is estimated and a 
joint sparse recovery step during which the active subcarriers 
are estimated. The OMP algorithm ll20ll . 1211 is tailored to 
recover the sparsity pattern of OFDM symbols. 

Paper organization: The rest of the paper is organized 
as follows. In Section |II] the system model is given and the 
problem is formulated. Sections |In] and |IV] are the main parts 
of the paper. In Section [Hll we first study the oracle-assisted 
case, then discuss the issues arising from very small or very 
large timing offset and introduce a solution for these cases. 
Then in Section |IV] we propose a joint compressive sensing 
method for identification of active subcarriers when the timing 
offset is known. We also extend the method to the case of 
unknown timing offset by employing a joint dictionary learning 
and sparse approximation algorithm to identify the active 
subcarriers. The performance of the proposed method will be 
studied by simulation experiments in Section |V] Conclusions 
are drawn in Section IVll 

Notations and Mathematical Preliminaries: Throughout 
this paper matrices and vectors are denoted by capital and 
small boldface letters, respectively. = denotes the equality 


and = denotes the definition. I and 0 represent, respectively, 
identity matrix and all-zero matrix. For a matrix A, [A]ij 
denotes the (i,j)-th entry of the matrix and As denotes the 
sub-matrix of A composed of only columns determined by 
index set S. The rank, column space, and null space (kernel) of 
A are denoted by rank(A), C(A), and /C(A), respectively.The 
set of integers from pto q > p (including p and q) is denoted as 
{p : q}. For a vector x, the (sub-)vector consisting of entries 
determined by index set § is denoted as xs. We denote the 
empty set by 0 and for a set S we denote its cardinality by |S|. 

Consider an M x N matrix A, the spark of A, denoted 
by spark(A), is the smallest number n such that there exist 
a sub-set of n columns of A which are linearly dependent 
E2\ . II23I . The following inequality describes the relaltionship 
between rank and spark ll22l : 

spark(A) < rank(A) -t-1. (1) 

If M > iV and A is full-column rank (and therefore there 
is no subset of columns that are linearly dependent) then we 
use the convention spark(A) = N + 1. On the other hand, if 
M < N and every M x Ad subset of A is invertible, then A 
is called a full-spark matrix ll24l . i.e. spark(A) = M -f 1. 

II. System Model 

Consider a primary user transmitting OFDM signals with 
No subcarriers out of which only K ^ No are nonzero. A 
cognitive secondary user receives the transmitted signal of PU 
and wishes to detect empty subcarriers to start communicating 
over them. Furthermore suppose that because of relatively 
high cost and energy consumption of conventional Nyquist- 
based Analog-to-Digital Converters (ADC), the cognitive user 
samples the OFDM signal at a sub-Nyquist rate to collect com¬ 
pressive samples using an Analog-to-Information Converter 
(AIC) m, 03- Matrix-wise, this can be described as 


last N — d columns. Let denote 1 = A^ iFd ^2 and Bd ,2 — 
Ad_ 2 I'd, 1 - The set of Ni, equations (|4|i can be rewritten as the 
following comprehensive form 

z = BdS -f w, (5) 


where 
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Bd ,1 
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0 ■ 
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0 
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III. Effects of timing offset 

Before proceeding with the proposed solution, we first 
investigate how the time delay in OFDM signal may affect the 
subcarrier identification process. We first remark that if there 
was no timing offset, i.e., d = 0, then we could have estimated 
the joint support of vectors s[n], n = l,...,Nt, by using 
the joint sparse recovery methods from multiple measurement 
vectors such as ||251 . 1261 . This is clear from (|4|i as when d = 0 
the first term in the right-hand side of the equation disappears. 
But when d > 0 we can no longer employ these algorithms 
as in this case each s[n], n = 1,..., Nb — 1 appears in two 
measurement vectors z[n — 1] and z[n] resulting these two 
measurements becoming correlated. 

Now consider the oracle-assisted case when we know the 
true set of active subcarrier indices §. Taking only the columns 
of Bd and entries of s corresponding to § and discarding the 
rest, @ can be rewritten as 

z = Bd,sss+w, (7) 


z[n] = Ax[n], n = 1,... ,Nb, (2) 


where Nb is the number of collected OFDM symbols, x[n] 
is the n-th noisy OFDM signal with unknown timing offset 
d, z denotes the sub-Nyquist sampled signal, and A is the 
M X N measurement matrix (notice that K < M < N), 
where N = No + Acp and Acp denotes the Cyclic Prefix (CP) 
length. The n-th OFDM symbol x[n] can be written as 


x[n] 


Fd^2s[n - 1] 
Fd,is[n] 


-I-w[n], n = 1,..., Ab, 


(3) 


where s[n — 1] and s[n] are information vectors at time 
instants n — 1 and n which are both A-sparse with identical 
support set §. Moreover, F = [Fo,F]^ where F is the 
No X No Discrete Fourier Transform (DFT) matrix defined 
as [F]p_q = , Fq is the submatrix of 

F consisting of only the last Acp columns of it, Fd,i is the 
submatrix of F consisting of its first N — d rows and Fd ,2 is 
the submatrix of F consisting of its last d rows. From Q, the 
sub-Nyquist sampled signal at time instant n can be written as 


where 
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where Bd,i,s — Ad,iFd,i,s. Bd, 2 ,s — Ad, 2 Fd,i^s, and Fd. 2 ,s 
and Fd,i,s are sub matrices of, respectively, Fd ,2 and Fd,i 
with columns corresponding to set §. Furthermore, sg = 
[ss[0]^,ss[l]'^,... ,ss[Ab]^]^ where ss[n], n = 0,...,Ab 
is the entries of vector s[n] corresponding to set § and 
B^”J, n = 0,..., Ab, denotes the (n -f l)-th block-column 
of bd.s- From (|8]l it is clear that vector sg can be recovered 
if B(i g is full column rank. The following theorems state the 
necessary and sufficient conditions for recovering sg. 

Theorem 1. A necessary condition for B^^ g having full 
column rank is that K < d < N — K. 


z[n] = Ad,iFd, 2 s[n - 1] -f Ad, 2 Fd,is[n] -f Aw[n], 

n = l,...,Ab, (4) 

where A^ 1 is the sub matrix of A consisting of its first d 
columns and 2 is the sub matrix of A consisting of its 


Proof: If d < K, then it is easy to see that rank(Bd,i) £ 
rank{Ady) < d < K and therefore rank(Bd,i,g) < K which 
means B(i_i g does not have full column rank. Since B^; 1 g is 
the only nonzero block in the first block-column of B^ g (see 
then matrix B^ g does not have full column rank too. This 












concludes that in order for Bd,s to have full column rank we 
need to have d > K. 

Similarly, if N — d < K, then it is easy to see that 
rank(Bd, 2 ) < rank{A.d, 2 ) < N — d < K and therefore 
rankled, 2 ,s) < K which mean that this time B,; 2 ,S does not 
have full column rank. Since B^; 2 ,S is the only nonzero block in 
the last block-column o/B^ §, then matrix B^ g does not have 
full column rank too. This concludes that in order for B^; § to 
have full column rank we need to also have d < N — K. □ 

Theorem 2. Assume that M > {Nb + 1)K/Nb and entries of 
matrix A are independent and continuous random variables. 
If K < d < N — K then B^ g will have full column rank with 
probability one. 

Proof: The proof is deferred to Appendix. □ 

The above theorems can be also explained with the fol¬ 
lowing intuition: when d < K the number of independent 
linear combinations of K unknown nonzero elements of vector 
s[0] appearing in the linear set of equation (|7]i, through linear 
mixture matrix Bd.i.s, is d < K, which makes it impossible 
to estimate these elements. Similarly, when N — d < K, the 
number of independent linear combinations of K unknown 
nonzero elements of vector s[Nb] appearing in the linear 
set of equation (|7]l, through linear mixture matrix Bd, 2 ,s, is 
N — d < K, which again makes it impossible to estimate 
these elements. For the rest of the vectors {s[n]}^T7^ this is 
not however the case as each of them appears twice in (|7]): 
once through B^ i_g and once through B^ 2 ,s- 

To solve the issue, when d < K, we cancel out the effect 
of unresolvable vector s[0] by projecting z[n], n = 1,..., Nb 
to the (M — (i)-dimensional null space of B^i i, denoted by 
to get 

zW = Pb,,iZW 

= PBd,iBd. 2 s[n]+P^_^^Aw[n], n=l,...,Nb. 

(9) 

The projection in (|9]l causes us to lose the information lying 
on the d-dimensional subspace spanned by the columns of 
matrix Bd,i. This is not however a big loss as in this case 
d is relatively a very small number. 

Similarly, when N — d < K, we cancel out the effect of 
unresolvable vector s[Ab] by projecting z[n], n = 1,..., Nb 
to the (M — A + d)-dimensional null space of Bd, 2 , denoted 
'’y Pb^.s’ to get 

= ^Bd,is[n - 1] + Pb^ 2 Aw[n], n = 1,..., Af,. 

( 10 ) 

This time, the projection in (fTOl i causes us to lose the informa¬ 
tion lying on the (A — d)-dimensional subspace spanned by 
the columns of matrix B^^ 2 - Again, this is not a big loss as in 
this case A — d is relatively a very small number. 

IV. Active Subcarrier Identification 
A. Known timing offset 

When the timing offset d is known, active subcarriers can 
be identified by exploiting sparse reconstruction algorithms 


tailored to the model described in the previous two section. 
When K < d < N — K, we tailor the OMP algorithm 11211 
as follows to estimate the joint sparsity pattern of vectors 
s[n], n e 0,..., Nb'. 


Algorithm 1: Proposed OMP-based joint sparse recovery 
algorithm for active subcarrier identification when K < 

d<N-K _ 

Data: z[n], n = 1,..., Nb, d. A, K, set z[0] = 0 
Result: active subcarrier index set § 
initialization: A: = 0, § = 0 
while A: < A do 

A: ■(— A; -f 1, 

Nb 

j = arg min (^di 

..,Na} n=l y 

§ = S U {_)}, 

Form Bd g as in ([8]l, 
sg = (B* *sB,.g)-iB^sZ, 
z^z - Bd.sSg, 

z[n] ^ z((n — 1)M -f 1 : nM), n = 1,..., Nb, 


z[n — 1] 
zfnl 


The above algorithm re-confirms why the cases of d < K 
and d > N — K should be excluded, as B^gB^ g can be 


rank deficient in those cases. Above hdj is the j-th column 


of matrix 


B£i,2 

Bd,i 


On the other hand, for the two cases of d < A and 
d > A — A, we can find the active subcarriers by applying the 
Simultaneous Orthogonal Matching Pursuit (SOMP) algorithm 
proposed in 1221 to the set of Nb equations (|9]l and (fTOl . 
respectively. 


B. Unknown timing offset 

If the timing offset is unknown, we cannot solely resort to 
the algorithm introduced in the previous subsection to identify 
the active subcarriers. However, in this case active subcarriers 
and timing offset can be estimated by alternating between the 
following two steps: 

• Step 1: Fix d and estimate s by solving the following 
optimization problem: 

s = arg min ||z — B^sjl^, 

S 

s.t. |§|=A (11) 

This step of the algorithm can be solved using the 
batch OMP algorithm introduced in Section IIV-AI 

• Step 2: Fix s to the value estimated in Step 1 and 
estimate d as 

d = argmin ||z - Bds||f, (12) 

d 

The above solution is in fact a joint dictionary learning and 
sparse approximation problem ll28l - ll^ . In the first step, we 
fix the dictionary B^; and estimate the joint sparse vectors 
s[n], n G 0,..., Nb, and in the second step we fix the joint 
sparse vectors s[n],n G 0,... ,Nb to the values estimated in 
the previous step and find the best estimate for dictionary B^ 
from the dictionary set {Bq, Bi, ..., Bat-i}. 













Fig. I. Probability of time delay estimation for an OFDM signal with number 
of active subcairiers K G {1,3,5}, compression ratio M/N = 0.5, and 
number of OFDM symbols Nf, S {10, 20}. 



Fig. 3. Probability of time delay estimation for an OFDM signal with 
compression ratio M/TV € {0.25,0.5,1}, number of active subcanlers 
K = 3, and number of OFDM symbols Nf, S {10, 20}. 



Fig. 2. Probability of active subcarrier identification for an OFDM signal with 
number of active subcanlers K G {1, 3, 5}, compression ratio M/N = 0.5, 
and number of OFDM symbols Nf, G {10, 20}. 

V. Simulation Results and Analysis 

In this section we study the performance of the proposed 
algorithm for joint estimation of delay and active subcarriers 
via some simulation studies. For simplicity, we consider an 
OFDM system with IFFT size 32 which means No = 32. The 
cyclic prehx length is chosen as A/p = A/,/4 = 8 and the total 
OFDM symbol length is thus N = No + Nee = 40. Subcarrier 
symbols are assumed to be drawn from a 16-QAM constella¬ 
tion with unit energy. The elements of the measurement matrix 
A are drawn from A/"(0,1) and its columns are normalized to 
have unit norm. We emphasize that these chosen values are 
just examples for carrying out numerical simulations, and are 
not as such related to the fundamentals of the proposed method 
in any way. Number of OFDM symbols taken for subcarrier 
identification is Nf, — 10 or Nf, = 20. The number of 
random trials is 5000 to gather sufficient statistics for reliable 
evaluation of performance. The OFDM timing offset in each 
run is chosen uniformly randomly as an integer in the interval 

[0,Ar-i]. 

In the first experiment we study the impact of the number of 


active subcarriers on both timing delay estimation and subcar¬ 
rier identiheation. The compression ratio is set to M/N = 0.5. 
We compute the results for three different choices of the 
number of active subcarriers K S {1,3,5} (chosen uniformly 
randomly within No total subcarriers) and two choices of the 
number of taken OFDM symbols Nt G {10,20}. Figure [T] 
shows the performance of the proposed algorithm for timing 
offset estimation of OFDM symbols. Moreover, Figure |2] 
illustrates the performance of the proposed algorithm in terms 
of probability of correctly identifying active subcarriers. The 
probability of active subcarrier identification is computed as 
the ratio of average of active subcarriers identified correctly 
over the total number of active subcarriers, where the average 
is taken over all random trials. As we can see from Figures [T] 
and|2]as the number of active subcarriers increases the OFDM 
signal becomes less sparse and therefore the performance 
of the sparse recovery algorithm which is a main part of 
the alternating method proposed in Section IIV-BI decreases. 
Furthermore, as expected, the performance of the algorithm 
improves as the number of collected OFDM symbols increases. 

It is noteworthy to mention here that an error in estimation 
of the timing offset does not necessarily result in an error 
in identifying active subcarriers. This can be observed from 
comparing corresponding points in Figures [T] and |2] For 
example, at the signal-to-noise ratio of —10 dB for K = 1 and 
Nf, = 20, from Figure[T]it is seen that that the timing offset has 
been correctly estimated with a probability of about 0.25, while 
Figure |2] shows that the active subcarrier in this case has been 
correctly identified with the probability of almost 0.77. This 
is concisely because for small errors in estimation of d, there 
is still a big chance to detect the subcarriers correctly. This is, 
in turn, because for two timing offsets d and d' = d-\-€ where 
e is a very small nonzero integer, there is a high correlation 
between a certain column of and the same column of B^/. 
A mathematical analysis of this observation will be given in 
an extended version of this paper. 

In the second experiment, the impact of the compression 
ratio on the performance of the proposed algorithm is studied. 
The number of active subcarriers is set to A = 3 and we 
choose the compression ratio from set M/N = {0.25, 0.5,1}. 
























Fig. 4. Probability of active subcarrier identification for an OFDM signal 
with compression ratio M/N G {0.25, 0.5,1}, number of active subcarriers 
K = 3, and number of OFDM symbols Nf, G {10, 20}. 


Figures [3] and |4] illustrate the performance of the proposed 
method in terms of probability of timing offset estimation 
and active subcarrier identihcation, respectively. As it can be 
observed the performance of the proposed method enhances 
as the compression ratio increases. This is also expected as 
the performance of the sparse recovery algorithm which is a 
main part of the alternating method proposed in Section IIV-BI 
improves with the number of measurements M. 

VI. Conclusion 

In this paper we proposed a technique for identihcation of 
sparse active subcarriers of an OFDM signal from sub-Nyquist 
samples in the presence of timing offset. We hrst studied the 
effect of timing offset and provided the sufficient and necessary 
conditions for signal recovery in the oracle-assisted case. Then, 
based on that, by tailoring the Orthogonal Matching Pursuit 
(OMP) algorithm, proposed a joint sparse recovery algorithm 
for identifying active subcarriers when the timing offset is 
known. Then we showed that for the case when the timing 
offset is unknown, the active subcarriers can be identihed 
by alternating between a dictionary learning algorithm, which 
delivers an estimate for timing offset, and the joint compressive 
sensing algorithm proposed for known timing offset case. The 
method can be used for detecting spectral opportunities in 
cognitive radio systems. The performance of the proposed 
algorithm was studied via simulations. 
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Appendix; Proof of Theorem[2] 

Before proving Theorem |2] we provide the following two 
lemmas which are used in the sequel to prove the theorem; 


1) spark(Fd_2) = min(Vo, d) + 1, 

2) spark(Fd^i) = min(A^o, N — d) + 1. 


Proof: We only prove 1) here. The proof of 2) is similar. 

F 

If Nq < d then we have matrix F^ ,2 = [^h] where 

-pH A F{2Aro_d+i:Aro}. Since the Nq columns of are 
linearly independent then so are the Nq columns of F 4^2 ond 
therefore 

spark(F(i_2) = -Vq + 1, when d > Nq. (13) 


On the other hand if d < Nq then F^.i = 
[F{7 Vo -d+i:No}]^■ Ffoni this it is easy to verify that 

Fd.i = [F{i:d}]^A, (14) 

where A is a diagonal matrix whose (i, i)-th diagonal element 
is defined as [A]j^i = where ujnq — 

Since [F^i.^/j.]^ is a full spark matrix h24'i . therefore its 
product with the non-singular diagonal matrix A is also full 
spark i l22l? .' i.e. 

sparkjFd.i) = d -f 1, when d < Nq. (15) 

The proof of Lemma is then concluded from and (123. □ 


Lemma 2. Consider the Q x P full column rank tall matrix $ 
and the QxT random matrix (P whose entries are independent 
continuous random variables. Then with probability one we 
have 


dim 


c($)nc(’®') 


P+min(Q,r)-min(Q,P-fr). (16) 


Proof: From m Equation (2.17)] we have 


dim 




rank{^) -f rank{'^) — rank{[^, ’4']) 

P -f min((5, T) — rank([^, ’S']). 

(17) 


If the entries of ’S' are independent continuous random vari¬ 
ables, then with probability one its columns are linearly 
independent of themselves and also of the columns of $ as 
long as P T < Q and therefore rank([$, ’S']) = P -\- T. 
When P -\-T > Q on the other hand rank([$,S']) = Q. 
Putting these in ( 1771 ) the proof is concluded. □ 

Now we proceed to prove Theorem |2l Eirst notice that 
matrix Bd.s is of size MN^ x K(Ni, -f 1) which makes it a 
tall matrix since M > (Nt, -\- 1)K/Ni,. Hence it is enough to 
prove that its columns are linearly independent. To this end, 
we hrst prove that if K < d < N — K then both 1 § and 
Bd, 2 ,s in ® are full-column rank with probability one. 

We give the proof for B^ 1 g here. The proof for B^ 2 ,s is 
similar. Erom ff2\ Chapter 4] we know that 

rank(Ad,iFd,2.s) = rankjA^^) - dim[C(Aji) n/C(FJ2 §)]• 

’ ( 18 ) 

Since entries of A are continuous random variables we have 


rank(Ad,i) = min(M, d), (19) 


Lemma 1. We have 













with probability one. Furthermore, from Lemma [T] we con¬ 
clude that 2 ,s has full-column rank since K < d. Therefore 
we will have dim(/C(Fj 2 §)) = d — dim(C(Fd, 2 ,s)) = 
d — rank(F(i 2 ,s) = d — K. From these and from Lemma |2] 
we will have 

dim[C(Aj^) n /C(FJ 2 ,s)] = d — K + min(M, d) — 

min((i, d + M — K) 

Since K < M (because of the sparsity condition) we will have 
min(d, d + M — K) = d and therefore 

dim[C(Aj;i) n /C(FJ 2 ,s)] = min(M, d) - K. (20) 

Inserting ( fT^ and (l20l i in (fTSl) we will have 

rank(Bd,i,s) = rank(Ad,iFd, 2 ,s) = K. (21) 

In a similar fashion we can prove that 

rank(Bd. 2 .s) = rank(Ad, 2 Fd,i,s) = K. (22) 

Now, from m Equation (8.3)] it is easy to see that 
rank([B^°g B^^gj) = 2K. Adding the block-columns of B^ § 
in ® one at a time and employing ifSTl Equation (8.3)] yields 


rank(Bd.s) = rank([B^°^ . 

■ • -“d.S J) 

= {Nb + 1)K 

(23) 

which concludes the proof. 

□ 
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